Spacetime Encodings IV - The Relationship between Weyl Curvature and Killing 
Tensors in Stationary Axisymmetric Vacuum Spacetimes. 
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The problem of obtaining an explicit representation for the fourth invariant of geodesic motion 
(generalized Carter constant) of an arbitrary stationary axisymmetric vacuum spacetime generated 
from an Ernst Potential is considered. The coupling between the non-local curvature content of the 
spacetime as encoded in the Weyl tensor, and the existence of a Killing tensor is explored and a 
constructive, algebraic test for a fourth order Killing tensor suggested. The approach used exploits 
the variables defined for the Backland transformations to clarify the relationship between Weyl 
curvature, constants of geodesic motion, expressed as Killing tensors, and the solution generation 
techniques. A new symmetric non-covariant formulation of the Killing equations is given. This 
formulation transforms the problem of looking for fourth-order Killing tensors in 4D into one of 
looking for four interlocking two-manifolds admitting fourth-order Killing tensors in 2D. 

PACS numbers: 



I. INTRODUCTION 

Very little is understood about the implications that 
the curvature of a spacetime manifold has for particle 
motion within the spacetime. In the context of Extreme 
Mass Ratio Inspiral (EMRI) gravitational wave observa- 
tions careful knowledge of particle motion around com- 
pact objects could lead to a spacetime mapping algo- 
rithm [l], 0] • This paper provides a framework in which 
the relationship between orbital invariants and the cur- 
vature expressed by the Weyl tensor can be explored. In 
particular, it formulates a constructive algebraic test to 
see whether a particular stationary axisymmetric vacuum 
(SAV) spacetime admits an additional invariant, a gener- 
alized Carter constant, assumed to result from a Killing 
tensor. Only SAV spacetimes that have two commuting 
Killing vectors dt and and thus can be generated from 
a complex Ernst potential are considered. 

The existence of a totally symmetric tensor Tr ai ,„ am \ 
of order to, which obeys the Killing equations 
T(ai—a m ;P) = implies that the quantity 



_ rp(ai—a m ) 



Pai ' ' ' Vol T : 



(1) 



is constant along a geodesic, and thus provides a con- 
stant of motion. In Eq. J!]). p a indicates the particle 
momentum. 

It was shown in Q that the condition that a SAV space- 
time admits a second-order Killing tensor places direct 
restrictions on the components of the Weyl tensor. In 
particular, it limits the Petrov type to D (HQ. The ap- 
proach to the problem of checking whether a particular 
SAV spacetime admits a Killing tensor of rank m is con- 
ceptually simple: all that is required is to formulate the 
condition that each component of the Killing tensor ex- 
ists. These integrability conditions result in a number of 
conditions on the Weyl tensor. Numerical experiments 
@) Q suggest that a large number of SAV spacetimes 
may possess a fourth-order Killing tensor. Now, a to- 
tally symmetric tensor of order four in four dimensions 



has 35 possible independent components and the Killing 
equations impose 56 conditions on the gradients of these 
components. In the direct check suggested, the Killing 
equations have to be satisfied in conjunction with the 
10 vacuum field equations. Furthermore, by writing out 
integrability conditions, many more equations and un- 
known fields are generated. The magnitude of the calcu- 
lation may make the notion of the practical implementa- 
tion of this idea seem absurd. Possibly for this reason, 
no literature on and no examples of fourth order Killing 
tensors in general relativity seems to exist 0, Q. This 
paper demonstrates how, with some finesse, it is possible 
to check whether SAV spacetimes admit a fourth-order 
Killing tensor, and construct its components. 

A number of ideas lead to the problem becoming 
tractable. Adopt the point of view that a set of con- 
stants of geodesic motion are equivalent to a coordinate 
system ideally suited to describing the motions of free 
falling particles in spacetime. Suppose now that an ob- 
server with his/her own inertial frame is conducting bal- 
listic tests and keeping track of the velocity and positions 
of a series of projectiles, in an attempt to discover exper- 
imentally what this special coordinate system is. Sup- 
pose further that the observer planned to one day com- 
municate the results of the experiment to someone else 
traveling through the same SAV spacetime, in an unam- 
biguous manner. He is concerned about how to orient 
the coordinate system of the experiment to do so with 
the greatest ease. A mathematician may suggest that 
both observer and the traveler orientate their coordinate 
systems according to the geometry of the spacetime, for 
instance along the principle null directions of the Weyl 
tensor, which in most spacetimes will be unique 0], or 
by selecting a transverse frame. Whether this advice is 
experimentally feasible is irrelevant for the rest of the 
paper, but a well chosen tetrad is introduced in Sec. ITTll 
and all quantities expressed with reference to it. 

The investigation [l], 0, [I[ into understanding the rela- 
tionship between curvature and geodesies was begun to 
explore the possibility of exploiting the algebraic prop- 
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erties of the solution generation techniques to provide 
a method of mapping spacetimes and cataloging their 
properties observationally. One of the simplest methods 
of mapping one spacetime onto the next, and one that is 
directly related to the underlying SL(2,R) symmetry of 
the SAV field equations, is the discovery of a Backland 
transformation. Developed by Harrison [9( and Neuge- 
bauer [lOj], the Backland transformations (BT) are used 
to generate a new solution from an existing SAV solution 
using a function known as a pseudo-potential (or in re- 
cent mathematical literature, an "integral extension") to 
carry out the mapping. The choice of variables used in 
these papers for the curvature quantities makes the field 
equations particularly transparent and easy to program 
using algebraic computer systems such as Mathematica. 
These variables are adopted in this paper. Ultimately, we 
can ask the question of what additional conditions, if any, 
should be imposed on the pseudo-potential so that for a 
given BT both the curvature components and the Killing 
tensor components are mapped to the next solution. 

It can be observed that the Killing equations for SAV 
spacetimes decouple in such a manner that a subgroup 
of the equations is equivalent to the Killing equations for 
a two manifold in the (p, z) plane. This allows us to use 
the geometric picture, derived from the theory of dynam- 
ical systems, of what a Killing tensor on a two-manifold 
actually represents 0], and to exploit the accompany- 
ing symmetries and analytic structure. The remaining 
Killing equations couple onto the two-dimensional decou- 
pled system in a "tree-like structure" . Found by inspec- 
tion, this property allows a large number of Killing tensor 
components to be sequentially eliminated by well-chosen 
integrability conditions. The remaining 10 components 
and their undetermined derivatives are arranged so that 
repeated differentiation introduces a minimal number of 
new variables, so the number of integrability conditions 
grows much more rapidly than the new variables. Finally 
the process terminates leading to a spectacularly overde- 
tcrmincd linear system for the components of the Killing 
tensor and certain of their higher order derivatives, con- 
sidered to be independent potentials. 

The coefficients of the potentials in the linear system 
are polynomials in the field variables and their deriva- 
tives (the highest derivative of the metric functions that 
is required is the fifth). Inverting the linear system and 
writing down the conditions that a solution can be found 
result in the Killing tensor components, given explic- 
itly as rational functions of the field variables and their 
derivatives, up to a scaling factor. The consistency con- 
ditions for the linear system determine the conditions on 
the field variables that are required in order for a Killing 
tensor to exist. The direct inversion of the overdeter- 
mined system will be termed the brute force method and 
is discussed in Sec. EU I argue that it is computation- 
ally feasible for existing computers to find the analytic 
answer. The possible caveat being that once the condi- 
tion on the field variables is written down, it may not be 
easy to identify the physical implications of the result. In 



practice, the main application of the brute force method 
would be to construct the Killing tensor components for a 
given spacetime, once they are known to exist. It should 
also be noted that the analysis is local, so given a space- 
time and enough local derivatives at a point the process 
described in Sec. I VII would provide an immediate check 
of whether the Killing equations are consistent at that 
point. A number of simplifications and insightbuilding 
special cases are suggested in IVII A[ using the general 
framework of the brute force method. 

A more elegant tack than the brute force method is to 
ask, given a particular SAV spacetime solution that ad- 
mits a Killing tensor, which of the BT or other solution- 
generation techniques preserves this property. This ques- 
tion is formulated clearly mathematically but not an- 
swered in Sec. IVII Bl 

The main effort in the calculation is to arrange quan- 
tities so that it does not mushroom out of control. The 
actual arrangement can only be understood by reading 
through the following pages in detail. The emphasis is on 
being able to characterize the relationship between cur- 
vature, the solution-generation techniques and geodesic 
motion, so the role of the Weyl tensor components is 
made explicit wherever possible. 

In the derivation for second-order Killing tensors given 
in [3J , it was noted that the problem of seeking a second- 
order Killing tensor in a 4D spacetime was equivalent to 4 
interrelated problems of seeking second-order Killing ten- 
sors on 2D spacetimes. This allowed the techniques used 
in two-degree-of-freedom dynamical systems [llj to be 
applied in seeking a solution. This observation prompted 
the search for a similar formulation of the 4D, fourth- 
order Killing equations considered in this paper. An 
equivalent symmetric formulation of fourth-order Killing 
equations having a similar property was found and is pre- 
sented in Sec. IVIIII In this alternative formulation, the 
equations are written in the form of the four fourth or- 
der Killing equations of a two-manifold with additional 
interlocking conditions. This formulation is particularly 
concise, and appears to offer the possibilty of an explicit 
analytic solution to the problem. Sec. IIXI concludes the 
paper with a brief summary of the main results, a few 
wild speculations, and a number of more serious com- 
ments about the outlook of future work. 



II. SAV METRIC AND FIELD EQUATIONS 

Any SAV spacetime with two commuting Killing vec- 
tors, d t and d$ can be represented by means of the Lewis- 
Papapetrou metric, 

ds 2 = e~ 2 ^ [e 2 -'{dp 2 + dz 2 ) + R 2 d<p 2 } - e 2 ^(dt - ud<P) 2 . 

(2) 

The metric functions can be determined entirely by solu- 
tions of the Ernst equation for the complex potential £, 

3t(£) V 2 £ = V£ ■ V£, (3) 
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where V = d. 



pp 



p p 



■ d zz and V = (d p , d z ) . The real 



part of the Ernst potential is 



=2V _ 



$l(£). Line integrals 



of the £ potential determine the functions 7 and u>. The 
function R is any harmonic function obeying the equation 
R.zz + R.pp = Oj and represents a coordinate freedom in 
this formulation of the metric. One choice that is often 
used is to set R — p. 

Equivalent to the Ernst formulation of the SAV field 
equations is a formulation introduced by Harrison J9fl and 
Neugebauer [l(| ■ The variables introduced in [i| [Tc| are 
most suited to performing Backland transformations and 
turn out to be directly related to the Ricci rotation coef- 
ficients computed for the tetrad, introduced in Sec. IIIII 
It is the Harrison-Neugebauer notation (referred to as 
M variables) that is adopted for the rest of the paper. 
This notation makes explicit the relationship between the 
Backland transformations, the choice of tetrad (Sec. 
and ultimately the Killing equations (Sec. IVI| . and cor- 
responding constants of motion. Furthermore, the field 
Eqs. (O expressed in terms of these variables are very 
easy and efficient to program for symbolic manipulation 
(The variable names differ from the original ones used.) 

Introduce the complex variables £ = l/2(p + iz) and 
C = l/2(p — iz) and define 



drR 

Mi=&7, M 2 = -|-, M 3 = 



M* 



a c7 , M* 



R 

d c R 



M* 



£ + £ 
d c £ 
£ + £'' 



=, Ma = 



d^£ 

£ + £' 
d ( £ 
£ + £' 
(4) 



Note that the '*' operation is not simply complex conju- 
gation, in fact All = M*, M 2 = M|, but M 3 — M| and 

W a = m 3 *. 

In terms of the M variables, the field equations can be 
expressed as 



M 1X = 
M 2X = 



A I 



3,C 



AL 



4,C 



M 2)? = 



-i(Af 3 Af 4 *+M 4 M*), 

-M 2 M*, 
1 
2 
1 
2 

-Ml + 2(M 1 AI 2 - M3M4). 



(M 2 m; + M 3 M|) - M3M3* + M3M4 ■ , 
(M 2 M 4 * + M 4 M 2 *) + .\ /;.!/, - M4M4 ) , 



(5) 



The remaining five field equations result from the 
complex-conjugate expressions of Eqs. ([5|) , which are also 
required to hold. It should be noted the field equations 
determine only certain of the derivatives of the M vari- 
ables. In particular, M 1 ^, M 3 M 4 and their complex 
conjugates are left free and encode the nonlocal content 
of the spacetime curvature. These quantities change as 
one solution is mapped onto another one. How they enter 
in the Weyl scalars is shown in Sec. IIIII 

For reference sake, we give the derivatives of the re- 



maining metric components, 

2d^ip = M 3 + M 4 , d^uj = R(M 4 - M^e' 2 ^ '. (6) 



III. COMPUTATIONAL TETRAD AND WEYL 
CURVATURE 

In order to understand the interrelationship between 
Weyl curvature, Petrov Classification, and the Killing 
equations, it is necessary to introduce a null tetrad. With 
the particular tetrad constructed below it turns out that 
the metric functions can be completely eliminated, and 
the problem written in terms of the M variables and the 
components of the Killing tensor on the tetrad. The first 
step is to define a local incrtial frame or tetrad frame by 
the basis vectors 



E\ = —e^dt + ue^c 
E 3 = ei-*dp, 



E 2 = e'^Rdfj), 

E 4 = e<-^dz. (7) 



The corresponding contravariant basis vectors are 



we* 

2 = ~R 4 ~R ' 



E 1 = e-*dt, 

E 3 = e^d p , E 4 = e^d z . (8) 

Using this inertial frame, define the null tetrad to be 



k=^=(E 1+ E 2 ), 1 = ±={E X -E 2 ), 

m=^=(E 3 +iE 4 ) 



1 



m= -y=(E 3 - iEi) 



(9) 



The Weyl Tensor coefficients are 

*o = C a0lS k a m l3 k' y m s , 

* 2 = C a /3 y sk a m f, m y l s , 
* 4 = C a p^sm a l l3 rri 1 l\ 



(10) 



and when expressed on this tetrad in terms of the Ms 
they become (setting V = e 2l ~ 2 ^) 

*° = 2l7(- 2M i M 4 + M^M* + 2Mf + M 4 * f ), 
1 



* 2 = —{M*Ma + M 2 M\ - 2M 4 M 4 *), 



f2 



* 4 = — (-2M1M4 + M3M4 + 2M 4 



M 



(11) 



and = \I>3 = 0. Note how the derivatives of the M 
variables that are not determined by the field equations 
enter the Weyl scalars ^4 and \l/o- The SAV solutions 
are of Petrov type D if the Weyl scalars expressed on 
this tetrad obey the identity \& *4 = 9 *|. In Sec.|V]it is 
shown, as an example of the general formalism how the 
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existence of a second order Killing equation implies that 
this condition holds. 

The Killing equations treated in Sec. HVI require knowl- 
edge of the Ricci rotation coefficients associated with the 
tetrad (fc, I, m, to). (See [5[ for a definition of rotation 
coefficients). Given Eq. (TTZ)) below, these coefficients 
serve as dials in the calculation which can be adjusted 
to generate simpler example problems, as discussed in 
Sec. M 

In terms of the M variables, the nonzero Ricci rotation 
coefficients are (recall that the 24 rotation coefficients are 
antisymmetric in the first two indices j a bc = — Ibac) 



Mi - Ml 



7123 



7131 



7132 



7231 



7232 



7343 



2M% 



Mn* 



Mn* 



MX - 2MS 



Mt 



7124 



7141 



7142 



7241 



7242 



M 4 - M 3 

2V / 2V ' 

M 2 - 2M 3 



2V2V 

M 2 - 2M 4 



where Mb = d^lnV 
coefficients are 



7344 = 



2Mi - M 3 




(12) 



2V2V 
- M4. Vanishing 



7121 = 7122 = 7133 = 7134 = 7143 = 7144 = 0, 

7233 = 7234 = 7243 = 7244 = 7341 = 7342 = 0. (13) 

These rotation coefficients clarify the relationship be- 
tween the M variables used in the BT and the Newman- 
Penrose tetrad formalism. Furthermore, if the spacetime 
is static then 7123 = 7124 = and if the gauge is chosen 
to be R = p, 7132 = 7231 = 7142 = 7231- 

Finally, before moving onto the Killing equations, we 
note that the constant matrix 77 that is used for raising 
and lowering indices on the null basis has four nonzero 



components rj 



12 



V 



21 



T, and rj 



34 



43 



f. 



IV. KILLING EQUATIONS 

On the tetrad basis, the Killing equations for a totally 
symmetric Killing tensor T of order to, or equivalently for 
the vanishing of the totally symmetrized intrinsic deriva- 



tive T, 



("1 



0, can be expressed in terms of the di- 



rectional derivatives and rotation coefficients as 



T, 



'■(I 



(ai"-a m ,6) Tfl Tj 7c(oi02 03"-Om^ 



(14) 



In the case of SAV spacetimes, the absence of functional 
dependence on the two coordinates t and </>, and the 
choice of tetrad, Eq. imply that the directional 

derivatives in the 1,2 (or k, I) directions vanish. As a 



result only directional derivatives in the 3, 4 (or to, to) 
directions need to be considered. For the tetrad cho- 
sen in Eq. (|9|), these derivatives can be expressed as 
to" = y=<9c and to" = ^=c^r. 

Both the directional derivatives and the rotation co- 
efficients contain a factor 1 / \ / 2V. When programming 
the check for the integrability conditions for the Killing 
components it is convenient to remove this last explicit 
dependence on the metric functions V — e 27 ~ 2 ^. To do 
so we multiply Eq. (fT4"|) by \J2V and let T denote a vector 
in which all the independent components of the Killing 
tensor T have been arranged, and we also let M and M* 
be four dimensional vectors M = [Mi, M 2 , M3, M4] and 
M* = [Ml , M 2 * , M| , M|] . 

Then the n Killing equations for a Killing tensor of 
order to can be represented as the system of equations 



C\T <: + C B T I = MC i D T 



M*C E T, 



(15) 



where i = 1, • • • ,n and the matrices C\, C B , C l D and 
C E , contain only integer entries. This equation makes 
explicit the coupling between the components of the 
Killing tensor and the curvature of the spacetime, ex- 
pressed in terms of M variables. 

The C % matrices contain a lot of structure inherited 
from the Killing equations. We shall exploit this struc- 
ture to limit the number of computations that have to be 
performed to check whether a specific spacetime admits 
a Killing tensor of order to. In particular, only the cases 
to = 2 and to — 4 have thus far been considered. 

It can always be shown that the Killing equations for 
the components 27^. ../.n with fcj = {3,4} decouple from 
the larger system, and are the m'th order Killing equa- 
tions for a two-manifold with conformal factor V. This 
observation allows us to use the geometric understanding 
obtained from the field of dynamical systems about what 
these Killing tensors represent to identify appropriate 
variables and possibly useful coordinate transformations. 

The remaining components that contain the indices 
i = {1,2} are ultimately coupled to this two- manifold. In 
the to = 2 case, by writing out the appropriate integra- 
bility conditions [3] they too can be written in a form of a 
second-order Killing equation on a two manifold distinct 
from V . As is shown in Sec. I Villi this property persists 
for the to = 4 case, where four two-manifolds admitting 
a fourth-order Killing tensor arc in effect sought. 

If a Killing tensor of even order is considered, it can 
be shown by explicitly writing out the Killing equations 
on this tetrad that the equations for Killing tensor com- 
ponents such as TAABk that admit an odd number of 
h L = {3,4} and A,B = {1,2} indices decouple from the 
rest, and can be set to zero, without loss of generality. 
A more subtle argument using the Jacobi metric 0] and 
the symmetries of the additional invariant on the two- 
manifold yields the same result. 

In the next two sections, the two examples for to = 2 
and to = 4 are used to illustrate the ideas described 
above more concretely. These examples are also used to 
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demonstrate the reductions, in some cases empirical, that 
simplify the problem to a computationally manageable 
size. 



V. SECOND ORDER KILLING TENSOR 

In the second order (to = 2) case, there are six non- 
zero components of the Killing tensor that have to be 
considered and ten Killing (n = 10) equations that limit 
the gradients of these components. A summary of the 
structure of the Killing tensor components and the gra- 
dients determined by the Killing equations is shown in 
Table HI To avoid any confusion in notation, and to il- 
lustrate the general method, the Killing equations are in 
some cases written out in detail for this example. We 
begin with the second order version of Eq. (|14|) . 

V2VT {abJ) = 2 V nm V2V ln{ab T f)m . (16) 

The 10 Killing Eqs. (fl6|) can be divided into two groups. 
Four equations for the three Killing tensor components 
T(fe 1/ ! C2 ), with ki = {3,4} decouple from the rest and 
are the Killing equations for a two-manifold with con- 
formal factor V. After setting M B = d^(lnV) and 
M B = (In V) these equations are 

733,c = M%T S 3, T 34) f = --(T 33 £ + M B T 33 ), 

^44,C = MbTu, T M £ = --(T 44 . ? + M B T 4 4). (17) 

The two equations on the left of (fT7|) state that T 33 /V is 
an analytic function of £, and T44/V is an analytic func- 
tion of £ @. It is in principle always possible to choose 
a gauge or i?-function in @ such that two functions are 
equal to the identity, as was done in [3j. The equations 
on the right of (fTT)) can be viewed as defining the gradi- 
ent of the function T34. However, the gradient can only 
be considered to be valid if the cross-derivatives cancel, 
or equivalently if the integrability condition 

T 33,« + ( M BT 33 ) rc = T44.cc + (M B T 4i ) x (18) 

holds. If one differentiates Eq. (Tj"5)) with respect to £ 
and substitute in the field Eqs. ([5]) and Killing Eqs. (fTT)) . 
then d^T 33 can expressed in terms of a linear function of 

( c£T 33 , 9^44), where p = {0, 1, 2} and it is considered 

to be determined by the integrability conditions. This is 
indicated in the fifth column of Table HJ A similar argu- 
ment holds for 9^44. The functions M B and M B and 
their derivatives d^M B and d^M' B that are not defined 

by the field Eqs. ([5]) always enter as polynomial factors 
in front of the T variables. 

The remaining six Killing equations define the gradi- 
ents of the functions Tii, T 22 and Ti 2 (the first row entry 
on Tabled. Setting T = [Tu, T 12 , T 22 , T 33 , T 34 , T 44 ]* and 



S = [Tu, T12, T 22 ]* these equations have the general form 
S' l ( = MC l D T + M*C l E T, 
Si = MG D T + M*C E T, (19) 

where C and C are constant matrices with integer en- 
tries. The integrability conditions of the function S are 
surprisingly simple. They are linear functions in only 
(9^T 33 , 9^44) where p = {0,1} with coefficients poly- 
nomial in the M-variables (at most quadratic) and the 
first derivatives of M^, (at most linear). No other 
components of the Killing tensor enter the integrability 
conditions. 

It can be shown [3j that appropriate linear combi- 
nations of the integrability conditions for the functions 
T{ab): where {A, B} take on the values {1, 2}, and of Eq. 
(fT8f result in three more equations similar to (fl~8|) . but 
for a metric with conformal factor distinct from V. This 
creates the geometric picture that one is looking for four 
two- manifolds that admit second-order Killing tensors. 

By writing out the integrability conditions (henceforth 
IC) for the Killing tensor components whose gradients 
are fully determined by the Killing equations (the first 
two row entries in Table U) and differentiating some of 
these IC equations to fix the higher order undetermined 
derivatives of the remaining Killing components (third 
and fourth entries in Table[lJ column 5), one can generate 
a system of linear equations for the remaining unknown 
functions represented in column 6 of Table [U Taking ad- 
ditional derivatives of the ICs (column 4) and substitut- 
ing in the field Eqs. ([5]) and Killing Eqs. (fT7|) as well as 
the expressions for the derivatives fixed by the ICs col- 
umn 5 further increases the number of equations while 
keeping the number of unknowns T 2 constant. In do- 
ing so we build up a large overdetermined linear system. 
The number of unknowns, and thus an indication of the 
size of the linear system is given in column 7 of Table |U 
The coefficients in this system are polynomials in the M 
variables and their derivatives. The field equations and 
Killing equations, as written in terms of M variables and 
expressed on the tetrad presented in this paper, are very 
easy to program in Mathematica. 

The overdetermined system of equations so con- 
structed can be solved for the undetermined Killing ten- 
sor components (column 6) up to an overall scaling. (The 
tensor components that have been eliminated by writing 
out their ICs can ultimately be constructed from a line 
integral.) The consistency conditions that ensure that a 
solution can be found provide polynomial conditions on 
the M variables and their derivatives. 

In the case of the second-order tensor we find, the re- 
quirement, amongst others, that 

= * *4 - 9^ (20) 

or that the SAV spacetime is Petrov type D. The other 
requirements, namely separability in a certain coordinate 
system, and the explicit solution of the problem, are dis- 
cussed more fully in [3j. 
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Properties of Second Order Killing Tensors in SAV Spacetimes 


m — 2, n = 10, K.E. involving 6 Tensor Components 


Column 1 


Column 2 


Column 3 


Column 4 


Column 5 


Column 6 


Column 7 


Tensor 


Number 

Ul VyUXJLJL LJ • 


Derivatives 
fixed by K.E. 


Number of ICs 

v J c Ilv'I dtcU 


Derivatives 
fivpH hv TP 


Resulting Unknown 

Pntoviti ale ~t~ t~\ nDtopininD 
JT ULcllllalS LU i-J v ' HJ 1 HHIlfe! # 


Number 
of Pot. 


T(AB) 


3 


d c ■ 




3 








T 3i 


1 


d< ■ 




1 








T 33 


1 


d ( ■ 






«?■ 


6fT 33 , p = {0,1,2} 


3 


T44 


1 










a^r 44 , p = {0,1,2} 


3 


Total 


6 Comp 


10 K.E. 


4 I.C. 






Li = 6 



TABLE I: Details of the Second Order Killing Equations (K.E.). Integrability conditions abbreviated as IC. Indices k and 
i take on the values 3,4 or equivalently indicate m,m components. Indices A and B take on the values 1,2 or indicate Z, k 
components 



The total number in the lower-righthand corner of the 
tables in many ways represents the largest possible lin- 
ear system that may be required if the coefficients of the 
undetermined potentials are to remain polynomial in the 
M variables, and one is working in a general gauge. A 
number of simplifications exist that decrease this number 
and will be mentioned in Sec IVII Al In the case of the 
second-order Killing equations, we can always work in a 
gauge in which the conditions on the M variables can im- 
mediately be written down (effectively £2 = 0). A more 
careful treatment of the second rank Killing equations 
is given in Q and the coordinate systems, in which the 
metric functions are separable, are found and classified 
there. 

The requirement that a second-order Killing tensor ex- 
ist on a SAV spacetime is very restrictive. Eq. (TT5|) 
makes explicit the coupling between the curvature of a 
SAV spacetime contained in the M variables on the right 
hand side and the components of the Killing tensor. One 
can take the point of view that the n equations of (|15p 
can be used as a representation of the M variables. With 
this perspective a rough counting argument immediately 
implies that it is impossible to represent the eight func- 
tionally independent M variables in terms of the six T 
variables. At the very least, a larger representation, with 
more free functions in the T vector, is required, to en- 
compass general spacetimes. It is also not known whether 
all possible allowed M functions consistent with the field 
equations can be written in the form (fl"5"|) . A practical 
method of checking is proposed in the next section for 
the fourth-order case. 



VI. FOURTH-ORDER KILLING TENSORS 

The orbital crossing pattern of the numerical integra- 
tion of the geodesic orbits @, [1] , the pole structure of the 
Ernst equation [l3| , the algebraic structure of the Weyl 
tensor [3] , and the group structure of the Backland trans- 
formations 0], all indicate that an exploration into the 



existence of a fourth-order Killing tensor on SAV space- 
times is warranted. Possible experimental applications 
l] provide urgency. 

The analysis is considerably more complicated and 
computationally expensive than the second-order case. It 
is also highly nonlinear, as opposed to the second-order 
case which can be linearized as in l 3] . In this section we 
use the technique outlined in Sec. [V] to make the problem 
of checking for and constructing a fourth-order Killing 
tensor on SAV spacetimes, vulnerable to a brute force 
attack. It is argued that it lies within the range of cur- 
rent computational capabilities to analytically compute 
the gradients of the Killing tensor components and the 
integrability conditions on the M variables representing 
the curvature for all SAV spacetimes. In some sense an 
answer is assured, although it may be ugly. In subsequent 
sections, refinements of this approach are suggested that 
may yield deeper insight into the field equations, and 
allow computations of insightful special cases with less 
computer power. 

On SAV spacetimes, the 56 Killing equations for the 35 
independent components of a general fourth-order Killing 
tensor decouple into two groups. As mentioned in Sec. lIVI 
the Killing equations for components such as Taabu 
and TAk 1 k 2 k 3 , with an odd number of fcj = {3,4} and 
A, B = {1,2} indices, form an entirely separate group 
from those with even pairs of indices. The equations 
that result from the odd group have more restrictive con- 
ditions on the curvature, and are trivially solved by set- 
ting them to zero. The more general case involving only 
the even group is considered. This reduces the number 
of tensor components to 19, and the Killing equations by 
half. Table HTl lists the derivative properties of the "even" 
Killing tensor components. 

The fourth-order Killing equations under consideration 
in this Section are written out in full in Appendix [XJ 
Three groups of equations can be identified. The de- 
coupled group governing the derivatives T(fc 1 fc 2 fc 3 fc 4 ), with 
fej = {3,4}, are the 6 fourth-order Killing equations for 
a two-manifold with conformal factor V, Eqs. (|A4|) and 
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TABLE II: Details of the Fourth Order Killing Equations (K.E.). Integrability conditions abbreviated as IC. Indices k and i 
take on the values 3, 4 or equivalently indicate m,m components. Indices A, B, C, D take on the values 1, 2 and indicate I, k 
components 



(jASjl . Two of these Eqs. ([AT]) can be combined to elimi- 
nate the T3344 component (Table HT1 third row). The sec- 
ond group consists often Eqs. (|A1[) defining the gradients 
of the five T^a^AsAz) components where A; = {1,2}. 
These components can be eliminated in favor of five in- 
tegrability conditions (Table ITTI first row, and Eq. (|B1[0 . 
Finally there is a buffer group of 12 Eqs. ([A"2]) . (|A"3j) in- 
volving the gradients of the mixed terms TA 1 A 2 k 1 k 2 , six 
of which (namely Eqs. (1A2|) ) can be eliminated in favor 
of the the three integrability conditions for the Ta 1 a 2 34 : 
terms (Table HH second row, and Eq. (|B2[) ) . In each 
case where an integrability condition is written down for 
a potential mentioned above, it is absent in the resulting 
IC. 

The "buffer" group does not have its counterpart in the 
second-order Killing equations. Now, a valid solution to 
the fourth-order Killing equations can be constructed by 
taking the exterior product of two-second order Killing 
tensors. This type of solution to the fourth-order Killing 
equations can be considered to be reducible. For irre- 
ducible fourth-order Killing tensors, the "buffer" intro- 
duces a functional freedom that will necessarily include 
a greater set of SAV spacetimes than those admitting 
second-order Killing tensors. The currently open ques- 
tion is how large this set is. 

The nine ICs in column 4 of TableHH (Eqs. JEE) (fB2|) 
and (|B3|) ) and the 34 unknown potentials in column 6 of 
Table Ull now form a closed system to differentiation (pro- 
vided the field Eqs. ([5]), Killing equations represented in 
column 3, Table[[T]or Eqs. (fA"3j) and (|A"5|) and the higher 
order derivatives represented by Column 5 are substi- 
tuted whenever necessary). In particular if you consider 
a system of equations generated by the nine ICs, their 
first derivatives, all their second derivatives with respect 



to ( and C, and the derivatives <9^IC and IC, a 
completely overdetermined system for the L4 = 34 po- 
tentials has been constructed. This system may be ex- 
cessive in that more derivatives of the ICs than what 
is necessary may have been taken. However it gives one 
freedom in choosing the pivots during the the elimination 
of the Killing components as one moves toward finding 
the conditions on the M variables, so that the answer is 
assured. The coefficients of the Killing components are 
polynomials in the M variables and their derivatives. At 
most fourth-order derivatives of the M variables enter 
these equations. 

The elimination process of the Killing components is 
at this stage ad hoc, guided mainly by inspection and 
symmetries in the tensor indices. It is essential to choose 
linear combinations of the equations such that the pivots, 
which become denominators in the subsequent manipu- 
lations, remain as simple as possible. Following this pro- 
cess it is possible to eliminate all but the first derivatives 
of the T3333, T4444, T3444, and T4333 components and 
the ten basic components themselves (the four remaining 
Tk 1 k 2 k 2 k 2 and the six remaining T Al A 2 kk components), 
before it is no longer apparent how to choose the pivots 
without introducing complicated denominators and the 
polynomials in the numerator have become so large that 
Mathematica has trouble with memory. Thus currently 
one is in principle one 14-by-14 full-matrix inversion away 
from an analytic expression for the fourth-order Killing 
tensor components, in terms of the M variables and their 
derivatives. The conditions that the M variables have to 
satisfy can be obtained by factoring the resulting consis- 
tency conditions after the matrix inversion is complete. 
The polynomials in M and their derivatives are very 
large. However both the matrix inversion and the subse- 
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quent factoring of integrability conditions lie within the 
range of current computational capacity. They have to be 
performed only once to get explicite expressions for the 
Killing tensor components for all SAV spacetimes. The 
only caveat is that once the final factoring step has been 
performed, the resulting conditions on the M variables 
may not be in an easily recognizable or compact form. 
In the following sections we discuss a series of examples 
that exploit the framework just given but seek ways to 
simplify the computations. 



VII. SIMPLIFICATION OF THE GENERAL 
FRAMEWORK 

A. General considerations 

Currently the greatest obstacle that hampers the fur- 
ther reduction of the linear system is the complexity of 
the polynomial coefficients that multiply the tensor com- 
ponents. Since the eight M variables and their deriva- 
tives up to the fourth enter the coefficients, any reduction 
in the number of M variables is welcome. One somewhat 
ineffective reduction is the choice of gauge R = p which 
implies that M 2 = M$. In static spacetimes, we find 
M 3 = M 4 , M* = M* and T nkk = T 22kk , T nn = T 2222 
and T1112 = Ti 222 by a symmetry argument. This re- 
duces the number of integrability conditions by three, 
and unknowns by six and makes the brute force approach 
more accessible, at the cost of generality. 

The analysis performed in the previous section is local, 
so the existence of a fourth-order Killing tensor can be 
disproved by showing that the consistency conditions fail 
at a particular point in space. The large overdetermined 
system generated in the previous section can thus be used 
to show that there is no consistent solution to the Killing 
equations given a particular spacetime. 

Another simplification that can be imposed is the 
choice of a gauge in which T4444/V 2 = T 3S33 /V 2 = 1, 
which I shall call the Killing gauge. This choice reduces 
the number of unknown potentials by ten, however it 
does nothing to reduce the complexity of the coefficients 
of the remaining components. Choosing this gauge re- 
moves the power of the formalism to check for the exis- 
tence of a Killing tensor given a metric, but it may be 
useful in understanding the structure of the equations. 
This was the case in second-order Killing equations [3j, 
where an analogous gauge choice was made. 

However the most powerful application of the tech- 
nique appears to be, given a metric that is integrable 
and admits a fourth-order (possibly reducible) invariant, 
and a parameterized family of metrics associated with it, 
to check which of these metrics retain the property of 
integrability. The next subsection suggests an approach 
that uses this idea, and that could in principle be used to 
check the integrability of all SAV metrics. This method 
also gives more insight into the choice of the M variables. 



Checking Backland tranformations, a tale of 
two manifolds 



This section suggests an elegant method that exploits 
the framework outlined in Sec IVII and couples it to our 
knowledge of the solution generation techniques, in or- 
der to ascertain which SAV spacetimes admit a fourth- 
order invariant. It was in large part the elegance of 
the solution generation techniques, and in particular the 
Backland transformations (BTs) that originally moti- 
vated the choice of the M variables. This choice makes 
programming the check proposed in the previous sections 
computationally possible. The approach given here can 
be viewed as an alternative to the brute force computa- 
tion. 

BTs map a valid solution of the SAV field equations 
onto a another with a different Weyl tensor. The map- 
ping is carried out with a so-called pseudo-potential, and 
for a given starting spacetime only 16 of these trans- 
formations exist [9(. The basic idea is to start with a 
metric such as Kerr, which is known to admit a fourth- 
order Killing Tensor, and to ask which of the 16 possible 
BTs [|[ maintain the property that the new spacetime 
also admits a fourth-order Killing tensor. The main ad- 
vantage of working off an existing solution is that the 
derivatives of the M functions that are not fixed by the 
field equations are known and thus the complexity of the 
polynomials entering the overdetermined linear problem 
is reduced. 

A more mathematical statement of this approach is 
the following. Consider two manifolds (A4,g) and (M., g) 
with metrics g and <?, that obey the SAV vacuum field 
equations. Suppose that M admits a fourth-order Killing 
tensor and that Ai and M. are related by a Backland 
transformation. Namely, N = BN and N* = B*N* , 
where N = (M 3 , M 4 , M 2 ) and B and B* are the 3 x 3 
transformation matrices. The entries of these transfor- 
mation matrices are very simple rational functions of the 
pseudo-potentials q and A. The 16 possibilities are col- 
lated in [9( , where however the notation may vary slightly 
from that used here. Furthermore M\ is fixed by the last 
equation of (|5|). 

The potential A is chosen to be A = 
and as a result 



(u-0/(u + o, 



d\ = - [A(A 2 - l)M%d( + A" 1 (A 2 



l)M 2 d(\ (21) 



The pseudo potential q lies within the prolongation struc- 
ture and obeys the Ricatti equation [9( 



dq = 
+ A- 1 



9(1 + q\)MS -(q + X)Mt + ~(1 - q 2 )\M 2 



q(q + A)M 4 - (1 + q\)M 3 + -(1 - q 2 )M 2 ) 



(22) 



The gradients of q and A are thus determined in terms of 
the M variables of the known solution. Furthermore, the 
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derivatives of the M variables that are not fixed by the 
field equations are also explicitely known. An overdeter- 
mined linear system of Killing equations such as that sug- 
gested in Table HI1 can then be built and used to check the 
conditions on q that are required to maintain fourth-order 
integrability. This approach has not been implemented in 
full, however the result would allow us to very accurately 
quantify the size of the subgroup of SAV spacetimes that 
admit a fourth-order Killing tensor and would lead to a 
deeper understanding about the relationship between the 
solution generation techniques and the geodesic structure 
of the resulting spacetimes. 

VIII. ALTERNATE FORMULATION OF 
FOURTH ORDER KILLING EQUATIONS 

The geometric origin of the variables used in the previ- 
ous sections is obvious by construction: namely, they con- 
form to the components of the Killing tensor expressed 
on an orthonormal tetrad basis adapted to the underlying 
geometry of the metric. In other words, the tetrad is cho- 
sen to respect the symmetry imposed by the Killing vec- 
tors of the SAV spacetime, and turns out to correspond 
to the transverse frame of that spacetime. The other 
spacetime dependent quantities that enter the equations 
are the rotation coefficients associated with the tetrad. 

In this section I present an alternate formulation of the 
Killing equations, which was found by inspection. While 
the system of equations presented here and those of Ap- 
pendix [X] are entirely equivalent, the former has symme- 
try properties that highlight the difference between static 
and stationary spacetimes, allowing easy simplification. 
The formulation presented here has the additional fea- 
ture that the spacetime-dependent quantities that enter 
the equations are functions of the metric variables and 
their derivatives rather than the M variables previously 
used. The relative simplicity of the resulting equations 
also allows insight into some of the features of the solu- 
tions that are sought. Looking at the structure of these 
equations led to an ansatz for the Killing tensor compo- 
nents on the equatorial plane of static spacetimes. This 
ansatz appears to be correct, as shown in [14]. The new 
variables, denoted by P<i : j>, are linear combinations of 
the Killing tensor components on the transverse tetrad 
and of functions entering the metric. While there is some 
vaguely systematic way of searching for the new variables, 



they are best motivated by the fact they have been shown 
to work to aid explicite calculation in practice. The orig- 
inal intent was to find a set of variables that cast the 
Killing equations in as symmetric as possible a form. 

Based the derivation of the second-order Killing ten- 
sors given in [3| , I suspected that one could find a formu- 
lation in which the Killing equations would take on the 
form of four "interlocking" fourth-order problems for a 
two manifold. The formulation given here is the result of 
that search. A clearly defined algorithm for finding these 
coordinates does not exist, except to say grope roughly in 
this direction and if it is the right way, the answer should 
have the form suggested in Q for fourth order Killing ten- 
sors on a two manifold, if not continue groping. In what 
follows I thus resort to the annoying approach of giving 
the anzatz for the transformation without fully being able 
to disclose how it was obtained, or what the quantities 
refer to, except to say that the resulting equations are 
easier to solve and conform to some intuitive picture I 
had of their existence before they were found. 

All the properties of the Killing equations previously 
discussed in Sees. IIVI and I VII are inherited by the new 
system. The naming convention of the new variables 
displays some of this structure. For the variable name 
P<i:j>j the index j gives some indication of the differ- 
entiability properties of the variable. The designation 
j = 0, for example, indicates that the Killing equa- 
tions fully determine its gradient, and that it can be 
removed from this system to yield an integrability con- 
dition. There are nine of these variables, P<i : o> with 
i £ {—4, ■ • • , — 1, 1, • ■ • ,5}. The indices j — 3 and j = 4 
indicate that the Killing equations fix derivatives with 
respect to £ and ( respectively. The index i labels differ- 
ent variables with a fixed derivative structure. The final 
linear system is built up of the ten variables -P<i : j> with 
i £ {1 • • • 5} and j £ {3, 4}, and of higher derivatives. 

Without further hesitation, make the following ansatz: 
define the set of variables that are already decoupled as 
the fourth-order Killing equation of a two-manifold as 

T3444 =P<i-A>e 2l ~ 2 ^ , T3334 =-P<i:3>e 27 ~ 2 '^ ', 
73333 =P<5:3>e 4 ~'~ 4 ' 4 ' , X4444 ^P<5A>e 4 ^' 4 ^ , 
T3344 =P<-1:0>. (23) 

Furthermore, set the mixed components Tabij equal to 



J1233 


1 

~12 


g 2 7 -4^ / 


-4e^P<i :3 > + 


e 4 ^ (-3c 2 P <2:3> 4 


- 6wP< 3:3 > - 


K6P< 4:3 >) 4 


-3P 2 P< 2:3 >), 


P1244 


1 

_ 12 


e 2 7 -4^ ( 


-4e 2 ^P <1:4 > + 


e 4 ^ (-3^ 2 P <2:4 > 4 


- 6cjP< 3: 4> - 


K6P <4 :4>) 4 


- 3P 2 P<2:4>) , 


P1133 




1 2 7 -4V 
4 


(e 4 ^ (^P <2:3 > 


— 2ojP <3:3> — 2P <4:3> ) -J 


- 2ite 2 * (toP <2:3> 


— P<3:3>) + P 2 P<2:3> 


Pll 44 




1 2 7 -4V 

4 


(e 4 ^ (^ 2 P <2: 4> 


- 2wP <3:4 > - 2P <4 :4>) ^ 


- 2Pe 2 ^ (WP<2:4> 


-P<3:4>)4P 2 P<2:4> 
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T 2 233 = " -e 2 ^ (e** (W 2 P<2:3> - 2^P <3 :3> - 2P<4:3>) - 2Pe^ {CJP <2 :3> ~ P<3:3>) + P' ' P<2:3>) 

T 2244 = - X -e 2 ^ (e 4 ^ (lu 2 P <2 . a> - 2 W P <3:4 > - 2P <4: 4>) - 2Re 2 ^ (loP <2A> - P <3 . A> ) + R 2 P <2A> ) 
Ti 2 34 (~4e 2 ^P < . 1:0 > + e 4 ^ (w (2P<_ 3:0 > - ^P<_ 2: o>) + 2P<_ 4 :0>) + P 2 P<-2:0>) , 



71134 ( W ( 2P <-3:0> - ^P<-2:0>) + 2P<_4:0>) + 2Pe 2 ^ (P<_ 3 :0> - ^P<-2:0>) + P' (-P<-2:0>)) , 



C2234 



= i e - 2 ^ (e 4 ^ (W (2P < _ 3: 0> - ^P<_2:0>) + 2P<_ 4 :0>) - 2Rf?+ (P<-3:0> - ^P<-2:0>) + P' (-P<-2:0>)) . (24) 



Finally, set the components of the form T(AiA 2 a 3 A4)j — {1) 2} equal to 

e- 4 ^P< 1:0 > (P + + 4e-^P <2 :o> (P + e^uf + 8P <3: o> (R + e^uf 

o 

e 4 ^ 



Pi 111 TTT 



1112 



P, 



1222 



9 

64 



^[e 2 *P<4:0> (R + e 2 * 



T 2222 — — 
64 



u) + e -r <5: o>j , 

4 



e~**P <1:0> (P - e 2 ^a/f + 4e- 2 ^P < 2:o> (e 2 ^ - P)° + 8P <3:0 > (P - e 2 V 



~ Jq [e 2 ^<4:0> (e 2 ^ - P) + e 4 ^P <5:0 >] , 

^e- 4 ^P <1:0 > (e 4 ^ 2 - P 2 ) + ^- 2 ^P<-2:0> - |e- 2 ^P <2:0> (P- 2e 2 V) 



(P + e 2 ^) 5 



! ( ?e 2 ^^P <3:0 > - ^P<-3:0> 



(i2 + e 2 ^) - ^P<. 4:0 > + |e 2 ^P <4:0 > (P + 2e 2 ^) + ^P <5:0 >, 



9 

64" 



^e- 4 *P 



<1:0 > (e 4 V - P 2 ) + ^- 2 V < _ 2:0 > + ^e~^P <2:0> (P + 2e 2 ^) 



(R - e 2 ^ 



^e 2 ^P <3:0> - ^P<. 3: o> ) (e 2 ^ - P) - ^P<_ 4: o> + ^e^P <4: o> (2e 2 ^ - P) + ^P <5: 



T lla2 =P<.1: > + ^ e ^P<l:0> (i? 2 - e^C 2 ) 2 + 



32 

ie^P^O + ^W-P<2:0> 



0>! 



(e^cu 2 - R 2 ) 



e 2 ^ W P < . 3:0 > - ^P<3:0> (R 2 - 3e 4 V) - e ^P<_4: > + ^e 4 ^P <4:0 > + ^P <5: o>. 
8 16 16 



(25) 



r 



2e 2 ~> 



3P 2 



e 2 7 f R 2 e -*l, _ u 2\ 
/4 = ^P^ ^ 



After this horror of an ansatz, the resulting Killing equa- 
tions are much more amicable and concise. Substituting 
the ansatz in Eqs. flSSJ), (J23JI, and ||S5J), into the Killing 

Eqs. ([AT]) . (TA"2l . (TAlfl) . (TM| and jM]) yields, after simpli- Note that in E q s - © ^<5:4> is an analytic func- 
ncation, the following equivalent system of Killing equa- tion of C and indicates a gauge freedom still present 
tions; the terms whose gradients are not fully defined in the metric. Without loss of generality one can set 
and that enter the linear system used to construct the P <5 : 4> = p <5:3> = 1, with raiiishing derivatives further 
integrability conditions are given below 

1 

4 
1 



^<z:3>,C = /,-/' 5:3 " ^P <5 ..3>X *€{!••• 4}, 



<i:4>,C 



f,J' 5:4 - ^flP<5A>X » G {1---4}, 



P<5:3>X = P <5:4>,C - - 



(26) 



The functions fi entering these equations are defined in 
terms of the metric functions as 



simplifying the expressions (J26J) to 

P<,:3>, C = -fi,o ~ - -he i G {1 • • • 4}. (28) 

The equations governing the field components whose 
gradients are fully described are now given by 

4 2 

P<-1:0>X = ~ g- P <l:3>/l,^ ~ g h P <l:3>X> 



fi = e 2 ^ = V, h 



2e 2 ^ 

3P 2 "' 



<-l:0>,C _ ~~ 3' P <l:4>/l,C - -^f\P<VA>Xi 



<-2:0>,C — - 2P<1:3>/ 2 ,C _ ^ P <l:i>X 
- 2/ 1 ^P<2:3> - /l- P < 2 :3>,C' 
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p 

<-2:0> X 


— — 2P<1 : 4> /2,C 


— /2-P<l:4>,( 




— 2 /l : £-P<2:4> 


— /l-P<2:4>,Ci 


P<-3:0>X 


= — 2i 3 <i : 3>/ 3 £ 


— y 3-f ^- 1 7 
v. 1 .0^ 




— 2 /j —P<3:3> 


^ A <vi:o>,C' 


P 

<-3:0>,C 


= — 2P<i : 4> /3,£ 






— 2 /l,^-f > <3:4> 


— /l-^<3:4>,(! 


-P<-4:0>,C 


= — 2i 3 <i : 3>/ 4 ^ 


<l.o>,C, 




— 2/^ ~P<4:3> 


r i P !*■ a . q 7" 1 

J A <4: o>,(, " 


P 

<-4:0>,t; 


= — 2-P<i;4> J4lC 


— J4-P<1:4> ) ^ 




— Z 1 1 f r<? A-A^ 


— M P 4- ■ 4- ^> ■> 


-P<1:0>,C 


= — 8P<2:3>f 2 X 


- 4/ 2 P <2:3> ,- 


P 


= — 8P<2:4>/2,C 


- 4/ 2 ^<2:4>,C 


P<2:0>,( 


' 1/' ±3 /,,- 


+ 2 h P <2:3>X 



+ 4/ 2!? P<3:3>+2/ 2 P <3:3>i? , 

P <2:0>,C = + 4P <2:4>/3,C + 2/ 3 P<2:4>,C 

+ 4/ 2 , C P< 3 :4> +2/ 2 P <3 :4>,C, 
-P<3:0>,C = + 2P<2:3>/ 4j 7 + /4-P< 2 :3> i 7 + 2/ 2 7-P<4:3> 

- 4P<3:3>/ 3; 7 ~ 2 f3 P <3:3>X + / 2P <4:3>,C' 
^<3:0>.C = + 2P< 2 :4>/4 j( ; + /4-P< 2 :4>,C + 2/ 2 ^P<4:4> 

- 4P<3:4>/3,C ~ 2/ 3 P <3: 4 > ^ + / 2 P<4:4>,C7 
f<4:0>,C = - 8P<3 : 3>/ 4 7 — 4/4P <3 . 3> 7 

~ 8 / 3 ,C P < 4:3 > ~ 4 /3-P <4:3>: 7, 
^<4:0>,C = _ 8P<3 : 4>/4 j( j — 4/4P <3: 4 >;( j 

- 8/3^P<4 : 4> - 4/ 3 P < 4 :4> ^, 

^<5:0>,c = — 8P<4 :3 >/ 4 7 — 4/ 4 P <4:3> 

P <5:0>,C = ~ 8P<4:4>/4 i( ; - AfiP <4:A> X . (29) 

An additional advantage of this formulation is that the 
terms on the right-hand side of the gradients do not con- 
tain any of the P<i : o> variables whose gradients are being 
defined. Thus, unlike the case treated in App. [AJ it is 
immediately obvious that the ICs contain only the ten 
variables P<i:3> and P<; : 4>. 

In the event that the spacetime is static rather than 
stationary /s = 0, and the following P<;.j> variables can 
be set to zero: 

-P<3:3> = -P<3:4> = -P< 2 :0> = -P<4:0> = P<-3:0> = 0. 

(30) 

In the general SAV case, the integrablilty conditions for 
this formulation can be generated by cross-derivatives of 
Eq. (f/T9"|) . Before continuing, however, it is useful to de- 
fine the following differential operators, which illuminates 
the structure in these equations. Let 

fP(i,j) = -'/,_/' , : ; + 3fi,cP<rA>,C + I 1 ' M 

TV(i,j) = 2/,_ + 3/ a P <J:3>j 7 + /,/',:;; 

(31) 



The I.C.'s for the fields {P<-i-.o> , P<i-.o> , P<5-.o>} can be 
written as 



TV(i,i) = TV(i,i) (32) 

where the index i takes on the values i E {1,2,4}. 
Note that these ICs are equivalent to the ICs on the 
conformal factor of a two-metric in a two-dimensional 
spacetime that admits a fourth-order invariant 0]. 
Thus in an abstract sense one is looking for at least 
three distinct interrelated two-metrics that admit fourth- 
order invariants in addition to satisfying the ICs gen- 
erated by the remaining fields. The ICs for the 
fields {P<- 2 :o>, J P<-3:0>,-P<-4:0>, J P<2:0>,-P<4:0>} take 
the form 



TV(i,j)+TV(j,i) = TV(i,j)+TV(j,i), (33) 

where the index pair takes on the values 

= {(1,2), (1,3), (1,4), (2, 3), (3, 4)}. The remaining 
IC for P<3:0> takes the form 

2TV(3, 3)-TV(2, 4) - TV(A, 2) = 

2.FP(3,3) - TV (2, A) - TV (A, 2), (34) 

which slightly breaks the symmetry of the previous two 
sets of expressions (|32p and (|3"3"|) . However, if in addition 
to the Killing equations, / 3 and P<3 : 4> obey the integra- 
bility conditions of a two-manifold with a fourth-order 
invariant, namely Eq. (|32|) with i = 3, the symmetry is 
beautifully restored. Eq. (f3"4"|) provides the missing pair 
— (2,4) in (j3"3"|) . The existence of a fourth-order 
invariant on a SAV spacetime then becomes synonmous 
with four interlocking two-manifolds admitting a fourth- 
order invariant (Eqs. (|3"2"|0 with additional conditions 
(Eqs. O)). 

These ICs coupled with Eqs. ([26]) , differentiated a few 
more times can also be used to build a large linear system 
to perform a brute force check. It does however turn out 
that they are of greater use as a framework for guessing 
the approximations of the invariant, whether or not it 
actually exists. A means of doing this systematically is 
explored in greater detail in [l5| and (l4| . 

IX. CONCLUSION 

In many ways the geodesies of a spacetime can be 
viewed as one of its most fundamental descriptions. The 
nature of the self-contained paths through space that 
massive particles favor gives us a direct observational 
characterization of the inertial field and of the mass dis- 
tribution within the spacetime. The field equations of 
SAV spacetimes are completely integrable. If the inte- 
grability properties are inherited by the geodesies of the 
manifold and that relationship can be made explicit; the 
practical application of a nonperturbative description of 
geodesic motion in SAV spacetimes can possibly be im- 
plemented in an experimental environment. 
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This paper provides a framework by which the relation- 
ship between constants of motion, expressed as higher- 
order Killing tensors, and the curvature content of the 
spacetime can be quantified. To do so most elegantly 
it is suggested that the solution generation techniques 
themselves be exploited to simplify the calculation. This 
will also increase our understanding of what these trans- 
formations imply physically for the spacetimes they are 
mapping. If the ability to maintain the existence of a 
fourth-order Killing tensor for certain of the Backland 
transformations of a given spacetime can be quantified, 
it will give some insight into the size of the subgroup of 
SAV spacetimes that admit a fourth order Killing tensor. 

The main sticking point at the moment is the the lack 
of a good method of eliminating all the Killing tensor 
components to obtain the conditions for the existence of 
the fourth-order Killing tensor on the M variables with- 
out building up excessively large polynomials. It is hoped 
that an elegant method of inverting the linear integrabil- 
ity matrix can be developed instead of the current ad hoc 
approach for choosing pivots. Failing that, however, the 
brute force method will eventually yield a result. 

If all SAV spacetimes admit a fourth-order Killing 
tensor, then there is a direct relationship between the 
geodesic structure contained in the Killing equations and 
the algebraic structure of the Weyl tensor. It further im- 
plies that the pole structure of the Ernst equation can 
be understood in terms of the poles of the analytic func- 
tions hidden in the Killing equations. The field of two- 
degree-of-freedom dynamical systems, with its paucity 
of examples of systems admitting a fourth-order invari- 
ant will then gain the solution generation techniques for 
Ernst's equations and the potential of a bi- infinite scries 
of examples generated by the HKX transformations and 
other solution generation techniques. Moser presented 
a geometric picture of how integrability on the Jacobi 
ellipsoids arises, and of the relationship to quadrics, by 
means of a very simple geometric construction. Neuge- 
bauer [10( pointed out that each axisymmetric stationary 
vacuum field corresponds to a minimal surface on a hy- 
perbolic paraboloid x + y + u 2 + v 2 — w 2 — embedded 
in a five-dimensional pseudo-Euclidean space. He further 
ventured that the Backland transformations were map- 
ping minimal surfaces into minimal surfaces. It is possi- 
ble that Moser's picture could be generalized to quartics 
which have their origin in the Lorentz boosts and that us- 
ing Neugebauer's observation a similar geometric picture 
could be built. This paragraph is speculative, it could list 
a lot of wishful thinking. In fact numerical evidence of 



chaotic behavior exists in some regions of phase space for 
the SAV spacetimes @ . This paper provides an algebraic 
test for the validity of these ideas. 

The alternative symmetric formulation of the Killing 
equations given in Sec lVIIII has several advantages, de- 
spite its currently obscure mathematical origin. By con- 
firming the hunch that the fourth-order Killing equations 
in 4D spacetime can be written in the form of four fourth- 
order Killing equations of a two-manifold with additional 
interlocking conditions, it gives greater insight into the 
structure of the Killing equations in 4D spacetimes. The 
high degree of symmetry of these equations and their 
compactness eases the way for further analytic explo- 
ration, increasing the possiblity for their analytic solu- 
tion, which will be attempted in [l5j . Furthermore, the 
simplifications that result from the additional assump- 
tion that the spacetime is static become immediately ap- 
parent, eliminating a large number of variables from the 
problem. This fact, coupled with an additional assump- 
tion of equatorial symmetry, allows one to obtain an ap- 
proximation for the Poincare map on the equatorial plane 
of these spacetimes. This program has been carried out 
for the Zipoy-Voorhees metric [HI and the results agree 
very well with numerical simulations of the orbits. Find- 
ing a more systematic derivation of the variables used 
in this formulation beyond just being the product of a 
fortuitous postulate should be an interesting exercise in 
itself. 

When gravitational-wave observatories such as LIGO 
and LISA mature and allow us to probe the nether re- 
gions of spacetime around compact objects, we stand to 
learn a great deal. It is hoped that this and the preceding 
papers in this series, [l|, |2, 3] will provide a mathematical 
framework in which to discuss practical algorithms for 
mapping spacetimes. So doing to facilitate optimally de- 
coding the information gleaned from gravitational wave 
observatories. 
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APPENDIX A: FOURTH ORDER KILLING EQUATIONS IN A TRANSVERSE TETRAD 

Killing tensor components confined to the directions set by the SAV Killing vectors have their gradients fully defined 
by the Killing equations. These equations, expressed in terms of the directional derivatives and rotation coefficients 
are given below, 

Tim, 3 = 2Tiin7i3i — 4Tni27i3i — 12Tn347i3i — 4Tim7i32 — I2T11337141 — 2T11117232, 
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Tim^ = —12Ti 1447131 + 27ii ii7i4i — 4Tin 2 7i4i — I2T11347141 — 4Ti 1117142 — 2T11117242, 

72222.3 = — 2T22227131 — 4T22227132 ~ 4Xi2227232 + 2T22227232 ~ 12722347232 _ 12722337242, 

72222.4 = — 2T22227141 ~ 4T22227142 _ 12722447232 — 4Ti 2 227242 + 2T22227242 ~ 12722347242, 

7 1 ni2,3 = (7'iii2 — 3Tn22 — 6T1234) 7131 — (4Tni2 + 6T1134) 7132 — 6T12337141 — 6T11337142 — (Tun + Tm 2 ) 7232, 
7 1 ni2,4 = — 6T12447131 — 6Tn447i32 + (T1112 — 3Tn22 — 6T1234) 7141 — (47\n2 + 6T1134) 7142 — (Tun + 7lh2) 7242, 

7l222,3 = — (7l222 + 72222) 7131 — (4Ti 2 22 + 6T2234) 7132 _ 6T22337142 + ( _ 3Tn22 + 7i 2 22 _ 6T1234) 7232 _ 6T12337242, 
7 1 1222,4 = — 6722447132 _ (Tl.222 + 72222) 7141 _ (4Ti222 + 6T2234) 7142 _ 6Ti2447232 + (~ 37u22 + 7i222 _ 6T1234) 7242, 
7 1 H22,3 = — 2(Ti222 + 72234)7131 — 4(Tii22 + 2Ti234)7l32 — 2T22337141 — 8Ti2337l42 — 2 (T1112 + 7ii34) 7232 — 2Tn337242, 
7 1 H22,4 = — 2T22447131 — 8Ti2447l32 — 2(Ti222 + 72234)7141 ~ 4(Tn22 + 2Ti234)7l42 — 2T11447232 — 2 (T1112 + T1134) 7242- 

(Al) 

The equations governing Killing tensor components with mixed indices TA 1 A 2 k 1 k 2 that have fully defined gradients 
are, 

Tll34,3 =2 ( — 71133,4 - 2Tii3 4 (27i32 + 7232 — 7l3l) ~ ^133(27142 + 7242 - 7141 + 27344)) 



(2Ti234 + 13344)7131 — (T1233 + 73334)7141, 
Tll34,4 =2 ( — 7ll44,3 - 2Tn34 ( — 7141 + 27142 + 7242) + T1144 (7131 — 27132 - 7232 + 27343)) 

— (T1244 + T 34 44)7i3i — (2T1234 + 13344)7141, 
T2234,3 =2 ( — T 2 233,4 - 2T 22 34(7l31 + 27132 — 7232) — T 22 33(7l41 + 27142 — 7242 + 27344)) 

— (2Ti234 + T 33 44)7232 — (T1233 + T 33 34)7242, 
T2234,4 =2 ( _ 7 1 2244,3 ~ 2 (+T2234 (7l41 + 27u 2 - 7242) +) - T 22 44 (7l31 + 27i 32 - 7232 - 27343)) 

— (Ti244 + T 3 444)7232 — (2Ti234 + T 33 44)7242, 

Tl234,3 =2 ( — Ti233,4 — 2T 2 2347l31 — T 2 2337l41 — 2T11347232 — T11337242 — 2T12337344) 
(2Ti234 + T 3 344)7i32 — (T1233 + 73334)7142, 



Tl234,4 ( _ 7i244, 3 — 722447131 — 2T 2 2347l41 _ T11447232 — 2T11347242 + 2T12447343) 

— (Ti244 + T 3 444)7i32 — (2T 12 34 + 73344)7142- (A2) 

Mixed components that have partially defind gradients are governed by the equations; 

2 

Tll33,3 = T1133 (7131 — 27132 — 7232 — 27343) — - (3T12337131 + T33347131 + T33337141) , 

2 

72233,3 = — g (3Ti2337232 + T33347232 + T33337242) — T2233 (7131 + 27132 — 7232 + 27343) , 
Tl233,3 = g ( _ 3T 2 2337l31 ~ 2T33347132 — 2T33337142 — 3T11337232 — 6T1233 (7132 + 7343)) , 

2 

Tll44,4 = — g (744447131 + (3Ti244 + T3444) 7141) + T1144 (7141 — 27142 — 7242 + 27344) , 

2 

T2244,4 = — - ((T44447232 + (3Ti244 + T3444) 7242)) — T 2 244 (7141 + 27142 — 7242 — 27344) , 

Ti244,4 = g (^2T 4 4447i32 — 3T 22 447l41 — 6T12447142 — 2T 3 4447 14 2 — 3T11447242 + 6T12447344) . (A3) 

The components orthogonal to the Killing directions T^k 1 k 2 k 3 k 4 ), h = {3,4} constitute a subgroup of equations that 
represent the fourth order Killing equations of a two manifold. Within this subgroup components with gradients fully 
defined are 

2 2 

73344,3 = — g (73334,4 + 2T33347344) , 73344,4 = — - (^3444,3 — 2T34447343) . (A4) 
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The remaining components with partially defined gradients are governed by the equations; 

73334,3 = — 2^3333,4 — T33337344 — 2T33347343, 23444,4 = — jT4444 > 3 + 744447343 + 2T34447344, 

73333,3 = — 4T3333 7343, ^4444,4 =47 4 44 4 7344. (A5) 



APPENDIX B: INTEGRABILITY CONDITIONS THE GENERAL FOURTH-ORDER CASE 



A large number of Killing Tensor components can be removed by writing down the integrability conditions (IC's) 
for all Killing tensor components whose gradients have been fully defined, namely equations (|A1|) , (|A2[) and (|A4|) . In 
particular the five tensor components Tim, 2iH2, T1122, T1222, T2222 can be completely removed from consideration 
using cross derivatives of Eqs. (|A1|) to yield the integrability conditions; 



(6T1244 + 2T 3 444)7j 31 + 3Tii44, 37131 + 2Tii447i3i,3 

+T1144 (-57i3i + (U7132 + 57232 - 47343) 7131) 

(6T1244 + 2T344 4 )7| 32 + 3T 2 244, 37232 + 2T 2 2447232,3 
+T2244 ((-57232 + 57i3i + U7132 - 47343)7232) 

' +3Ti244, 37131 + 3Tii44 ; 37i32 + 2Ti2447l31,3 
+3T 22 447l31 + 4T344471327131 + 9T114471317232 
+3T11447132 (47132 + 7232 ~ 27343 — 7131) 
\ +2Ti244 (-7131 + (107132 + 7232 - 27343) 7131) _ 

/ +3Ti244, 37232 + 3T2244,3 + 2Ti2447232, 37132 
+3Tn447f 32 + 4T344471327232 + 9T 2 2447l3l7232 
+3T2244 (7132 (7131 + 47132 — 7232 — 27343)) 

V +2T1244 ((—7232 + 7131 + IO7132 — 27343)7232) 

/ 3T2244, 37131 + 12Ti244, 37132 + 371144,37232 
+2(T 2 2447l3l,3 + ^11447232,3) 
+T3444 (87^32 + 471317232) 

+T2244 (7131 + (267132 - 7232 - 47343) 7131 ) 

+ 12Ti244 (47132 - 273437132 + 371317232) 
V +T1144 ((7232 — 7i3i + 267132 — 47343)7232) . 



(6T1233 + 2T3334)7^ 41 + 3T1133, 47141 + 2T11337141,. 
-T1133 (57i4i - ( 14 7l42 + 57242 1 



- 47344) 7141 

(6T1233 + 2T3334)7| 42 + 3T2233, 47242 + 2T22337242 

-T2233 ((57242 - 57141 - I47142 - 47344)7242) 

+ 3Ti233,47l41 + 3Tn33, 47142 + 2Ti2337l41,4 \ 
+ 3T 2 2337l41 + 4T333471427141 + 9T113371417242 
+ 3Tll33 (7142 (47142 + 7242 - 7141 + 27344)) 
-2Tl233 (7l41 - (107142 + 7242 + 27344) 7l4l) 

+ 3Ti233,47242 + 372233,47142 + 2Ti2337242,4 \ 
+ 3T11337I42 + 4T333471427242 + 9T 22 337l4l7242 
+ 3T 2 233 (7142 (7141 + 47142 — 7242 + 27344)) 
+ 2Ti233 ((-7242 + 7141 + 107142 + 27344)7242) / 

3T2233, 47141 + 12Ti233,47l42 + 3Tn33, 47242 ' 
+2(T22337l41,4 + 711337242,4) 

+4T3334 (2jf 42 + 71417242) 

+12233 (7l41 + (267142 - 7242 + 47344 ) 7l4l) 



-127 



(47i 2 42 



(Bl) 



> -— 1233 (47142 + 273447142 + 371417242; 
-7ll33 ((-7242 + 7141 - 267142 - 47344)7242) / 



Note that the five IC's of Eq. (|B1[) only contain first derivatives of the Killing tensor components. The three tensor 
components T1134, T2234, T1234 can be eliminated using cross derivatives of Eqs. (|A2j) to yield the integrability 
conditions; 



/ T3444 (-67^31 + 2 (127132 + 37232 - 77343) 7131 + 67131,3) \ 
■f 3T1144 (7f 31 + (-47132 + 27232 + 37343) 7131 + 67132 + 7232 + 27343 ~ 7131,3) 

■f3Tii44 (471327232 + 7232,3 — 871327343 — 372327343 — 27343,3) 
-6T1244 (7l31 - ( 4 7l32 + 7232 ~ 37343) 7131 - 7131,3) + 6T 22 447l31 
\ 1073444,37131 + Tii44, 3 ( — 67131 + 127i 32 + 67232 — 97 343 ) + 12Ti 2 44,37l31 + 37ii 4 4,3 3 / 

f 2T3444 ( — 37I32 + 371317232 + 1271327232 — 773437232 + 37232,3) \ 
+3T 22 44 (7l31 + ( 4 7l32 + 27232 - 37343) 7131 + 67132 + 7232 + 2 7343 + 7131,3) 
+3T 2 244 ( — 471327232 — 7232,3 — 871327343 + 372327343 — 27343,3) 
+6T1244 ( — 7232 + 71317232 + 47i 32 7232 — 37 3 4 3 7 2 32 + 7232,3) + 6T11447I32 
\ +10X3444,37232 + 12Ti244, 37232 + 3T 2 244,3 (27131 + 47132 — 27232 — 37343) + 37 2 244,33 / 



Interchange 
ind. 3 and 4 



Interchange 
ind. 3 and 4 
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' 2T 34 44 (97132 - 1073437132 + 971317232) ^ 
+ 3Tii44 (7332 — 71317232 + 471327232 — 373437232 + 7232,3) 

+6T1244 (37i32 ~ 473437132 + 7| 43 + 371317232 - 7343,3) 

+ 3T 22 44 (7l31 + ( 4 7l32 - 7232 - 37343) 7131 + 7131,3) 
\ +10T3444 i3 7i32 + 6Tii44 :3 7232 + 3Ti 2 44,3 (47132 ~ 37343) + 6T 2 244, 37131 + 7l 2 44,33 / 

The structure of the three ICs of Eq. (|B2[) differ slightly from that of Eqs. (|Bip in that they contain second derivatives 
of the Killing tensor components. Finally the component T3344 can be eliminated using cross derivatives of Eqs. (|A4|) 
to yield the IC; 

2T3444 (7343 — 7343,3) ~ 3T3444 7 37343 + T3444 i 33 = 2T3334 (7344 + 7344,4) + 3X3334^7344 + T'3334 i 44 (B3) 

Eqs. (|Bip . (IB2|) and (|B3|) constiture the nine basic integrability conditions on the remaining fields used to build up 
the linear system of equations for the Killing tensor compon ents and their derivatives. 



Interchange 
ind. 3 and 4 



(B2) 
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